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ABSTRACT 

We study the constraining power on primordial non-Gaussianity of future surveys of the large- 
scale structure of the Universe for both near-term surveys (such as the Dark Energy Survey 
- DES) as well as longer term projects such as Euclid and WFIRST. Specifically we perform 
a Fisher matrix analysis forecast for such surveys, using DES-like and Euclid-like configura- 
tions as examples, and take account of any expected photometric and spectroscopic data. We 
focus on two-point statistics and we consider three observables: the 3D galaxy power spec- 
trum in redshift space, the angular galaxy power spectrum, and the projected weak-lensing 
shear power spectrum. We study the effects of adding a few extra parameters to the basic 
ACDM set. We include the two standard parameters to model the current value for the dark 
energy equation of state and its time derivative, wq, w a , and we account for the possibility of 
primordial non-Gaussianity of the local, equilateral and orthogonal types, of parameter /nl 
and, optionally, of spectral index ny NL . We present forecasted constraints on these parameters 
using the different observational probes. We show that accounting for models that include pri- 
mordial non-Gaussianity does not degrade the constraint on the standard ACDM set nor on 
the dark-energy equation of state. By combining the weak lensing data and the information 
on projected galaxy clustering, consistently including all two-point functions and their covari- 
ance, we find forecasted marginalised errors <x(/ NL ) ~ 3, cr(nf NL ) ~ 0.12 from a Euclid-like 
survey for the local shape of primordial non-Gaussianity, while the orthogonal and equilateral 
constraints are weakened for the galaxy clustering case, due to the weaker scale-dependence 
of the bias. In the lensing case, the constraints remain instead similar in all configurations. 

Key words: Primordial non-Gaussianity — Large-scale structure of the Universe — Cosmo- 
logical parameters. 



1 INTRODUCTION 

The accurate measurement of the cosmic microwave background 
(CMB) anisotropics by the WMAP satellite l |Larson et al.|201 1\ , to- 
gether with the observation of Type la Supernovae and the deep and 
wide mapping of the large-scale structure (LSS) of the Universe, 
as done by the SDSS survey (Aihara et al.| 201 1 ), have led to the 
standard model of cosmology, whose parameters are now measured 
with unprecedented accuracy. However, the simplest model that fits 
all observations, ACDM, has some shortcomings and poorly un- 
derstood phases, in particular the late-time acceleration (or dark 
energy) (Frieman et al. 20081 and the early-time inflationary phase 
( |Linde||2008 1. For these reasons, it was agreed that new genera- 



tions of surveys should be built to improve upon the current data. 
The group known as Dark Energy Task Force (DETF) (Albrecht 
|et al.|20 06 1 has then labelled current finished and ongoing surveys 
as stages I and II, while future experiments belong to stages III 
or IV depending on their timescale and forecasted power, such as 
the Dark Energy Survey (DE S) <|Abbott et al.|2005) and the Euclid 
satellite (Laureijs et al. 2009. 2011 1 respectively. 
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These surveys will also shed new light on other basic ques- 
tions, such as the origin of the large-scale sU'ucture and the physics 
of the early universe. It is currently assumed that the observed inho- 
mogeneities in the matter density were seeded by quantum fluctua- 
tions at primordial times. These fluctuations were then stretched by 
inflation (or some alternative theory) on super-horizon scales, and 
acted as seeds first for the observed CMB anisotropies, and then 
for the matter density fluctuations. The simplest possible model for 
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the primordial fluctuations consists of a Gaussian random field with 
nearly scale-invariant power spectrum (Bardeen et al. 1986). How- 



ever, many models exist for the inflationary phase i Lyth & Lid-| 
|dle||2009) , many of which would alter the statistics of the fluctu- 
ations, in particular introducing a non-Gaussian component ( |Chen| 
|2010| l. It is therefore important to put observational constraints on 
the presence (and, once detected, on the amount) of primordial non- 
Gaussianity (PNG), in order to distinguish between models of the 
early universe. This is usually done in terms of the bispectrum am- 
plitude / NL for a set of template shapes for the bispectrum of the 
perturbations (Liguori et al. 2010|) (see the Appendix [A] for a more 
precise definition). 

The most well-established method to constrain PNG is to mea- 
sure the 3-point statistics of the CMB temperature anisotropies 
l |Komatsu|2010[ >. Current limits from WMAP7 for the local, equi- 
lateral and orthogonal types of PNG are -10 < /*£ < 74, -214 < 

< 266 and -410 < /°' L ' < 6 at 95% c. 1., respectively |k<> 
|matsu et al.|201 1) >. This will be greatly improved by the Planck mis- 
sion, for which the expected uncertainty is of the order 0"(/nl) - 5 
(Komatsu & Spergel 2001), limited by cosmic variance. A comple- 
mentary way is to constrain PNG from LSS studies. In this case, 
one can exploit multiple effects. First, the abundance of the most 
massive structures existing at any given cosmic time strongly de- 
pends on the level of PNG ( |Lucchin & Matarrese|1988| |Pillepich| 
|et al.p OlO) although the effect is degenerate with other cosmo- 
logical parameters (Pillepi ch et al.|2011) . The recent detection of 
a few high-redshift clusters gives some evidence in favour of non- 



Gaussian models ( Hoyle et al.|201 1 1, even though no consensus has 
been reached ( Mortonson et al.|201 l[|Ffoyle et al.|2012) also due 
to the low number statistics. The abundance of weak lensing peaks 
has also been proposed as a method to measure PNG (Mari an et al.| 
|201 1| |. Moreover, in full analogy with the CMB case, any 3-point 
statistic of the galaxy distribution is affected by PNG ( Sefu satti &| 
Komatsu 2007, Jeong & Komatsu 2009). However, the situation 
is complicated by the overlapping effects of primordial and late- 
time non-Gaussianity, which is driven by the non-linear growth of 
structure and galaxy biasing (Komatsu et al. 2009). These effects 
can be disentangled only by a thorough understanding of the non- 
linear regime. Similar considerations apply to modifications of the 
matter power spectrum on small scales induced by PNG ( |Taruya| 
|et al.|2 008; Pille pich et al.|20T0) . This is why the recent discov- 
ery that PNG also generates a strong scale-dependence of galaxy 
biasing on very large scales (Dalai et al. 2008) has attracted partic- 
ular interest in the literature. It has been shown that measurements 
of PNG based on this single feature are currently competitive with 
CMB studies ( |Slosar et al.|2008||Xia et al.|20lT) and will remain 
so in the future ( [Afshordi &~T olley 2008 ; Ca rbone et al.|2010J . 

In this paper we want to investigate further these techniques, 
applying the Fisher matrix formalism to determine to what accu- 
racy we expect future surveys will constrain PNG using two-point 
statistics only. We will compare the constraints from the large- and 
small-scale effects of the LSS with those expected by the CMB 
Planck mission, extending the results by|Joachimi & Bridle|(|2010|>; 
|Fedeli & Moscard ini (2010 ); |Carbone et al.| ( |2010| >; |Wang et al.| 
( 20 1 1 ; Namikaw a et al.H201 1\ . In particular, our study differs from 
the previous ones as: the forecasts on PNG by Fedeli & Moscardini 
( |2010} only considered weak-lensing data and did not include vari- 
ations on any other cosmological parameter but /nl; the results by 
|Carbone et al.| ( [2010) only considered the 3D galaxy power spec- 
trum and only the local type of PNG, while |Joachimi & Bridle] 
l |2^T0) ; [Wang et al.| ( |20To) did not consider PNG. 

The plan of the paper is as follows: we will review the effects 



of primordial non-Gaussianity on the LSS in Section [2] while, in 
Section[3] we will describe how we calculate the non-linear power 
spectra (for galaxy clustering and weak-lensing studies) in terms of 
the cosmological parameters. We will then describe the future sur- 
veys considered in Section]?] present our basic forecasts in Section 
[5] and expand on them in Section [6] We shall finally conclude in 
Section|7] 



2 THE EFFECTS OF PRIMORDIAL NON-GAUSSIANITY 
2.1 Definitions 

The simplest single-field, slow-roll model for inflation gives rise to 
a nearly Gaussian distribution of the curvature perturbations ( or 
the Bardeen potential O at primordial times corresponding to a red- 
shift z,. This changes however in most generalisations: many mod- 
els, and especially multi-field inflation, produce non-Gaussianities 
(see e.g. the recent review by |Byrnes & Choi|2010| >. There are dif- 
ferent possibilities for a departure from a purely Gaussian distribu- 
tion. The most general expression of a deviation from Gaussianity 
at quadratic level can be expressed by a non-local relationship be- 
tween the primordial Bardeen's potential O and a Gaussian auxil- 
iary potential ip. In real space ( |Schmidt"& Kamionkows ki|2010| l 



<5(x, z,) = <p(x, z.) + (/nl * W * (f * ip) (x, z.) , 



(1) 



where the asterisk denotes convolution and W(y, z) is a kernel 
whose form describes the type of non-Gaussianity considered. This 
is often called "the CMB definition", as opposed to "the LSS def- 
inition", as it is written at early times z». Here the / N l function 
quantifies the amount of PNG at first order, and it represents the 
first relevant deviation to measure, i.e. the skewness of the pertur- 
bations' probability distribution at a given length scale. Subsequent 
contributions at higher order, e.g. the kurtosis contribution £nl are 
also expected in many theories, but will not be considered in the 
following. In the simplest models, / N l does not depend on scale, 
and thus it simply assumes a constant value in Eq. Jl). In most of 
the following we will assume this simplification, but we will extend 
the analysis to the scale-dependent case in Section|63] 

We can define the power spectrum and bispectrum of the po- 
tential in the usual way as 



<0(k)0(k')) 
<®(k) d(k')0(k")> 



(27r) 3 <5 D (k + k')/\t,(£) 

(2/r) 3 6 D (k + k' + k") B^k, k', k") , (2) 



where the tilde denotes Fourier transformation. Neglecting sub- 
dominant corrections from the 4-point correlator (trispectrum), we 
can also assume P$,(k) = P<p(k) at leading order in / N l- Then it can 
be seen that applying the definition of Eq. <JTJ, the kernel W defines 
the relationship between power spectrum and bispectrum as 

B^(k, k',k") = 2/ NL [W(k, k') PtoQc) P*(k') + 2 perms.] , (3) 

assuming constant / N l- In the simplest case of PNG of the lo- 
cal type, W = 1 and the bispectrum peaks for squeezed triangles 
(k <k k' ~ k")\ besides this case, we will consider in the follow- 
ing also the equilateral configuration, for which the bispectrum is 
maximum for k ~ k' ~ k" , and the orthogonal type, which was 
constructed to be orthogonal to the previous two types. These con- 
figurations or a mixture thereof can be produced under different 
inflationary scenarios, see e.g. Babich et al. (2004); Senatore et al. 
l |2010} for more details. The bispectra for the three cases are given 
in the Appendix [A] 
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The analysis of CMB bispectra from the Planck satellite is cur- 
rently considered to be the most promising tool for distinguishing 
between these scenarios (Fergusson & Shellard 2009). In this pa- 
per we will consider the three configurations separately, although 
an analysis of a general linear combination of the three modes is in 
principle possible. 

The primordial Bardeen potential (equal to the gravitational 
potential with the opposite sign for sub-horizon modes), is then 
related by the Poisson equation to the total matter density perturba- 
tions 6. At linear level 



<J(k, z) = a{k, z) ®(k, z«) =s a(k, z) #(k, z.) , 
where 

2c 2 k 2 T(k)D(z) g(0) 



a(k,z) 



m„,H 2 g(z,) ' 



(4) 



(5) 



and we have introduced the linear growth function D(z) (nor- 
malised so that D(0) =1), the transfer function T(k), and the po- 
tential growth function g(z) oc (1 + z)D(z). We can therefore write 
the tree-level matter power spectrum 



P (k,z) = a 2 (k,z)P< b (k,z.) - a 2 (k,z)P v (k,z«). 



2.2 LSS and primordial non-Gaussianity 



(6) 



The large-scale structure of the Universe is commonly described in 
terms of different tracers. In high-density regions of the underlying 
dark-matter density contrast S, we observe the formation of dark- 
matter haloes, whose distribution is described by the field S/,. We 
expect that galaxy formation occurs within these haloes, thus pro- 
ducing a galaxy density field 6 g . On very large scales and to first ap- 
proximation, the linear bias coefficients b/, and b g relate the differ- 
ent density fields as 6i, — b/, 6 and S g - b g d. From the observational 
side, the dark-matter distribution 6 can be directly observed with 
weak-lensing studies, while S g can be measured by mapping the 
galaxy distribution. The most natural 2-point statistic that can be 
observed is then the power spectrum of these density fields which 
can be defined in analogy with Eq. l|2j. Therefore, the galaxy and 
matter power spectra approximately satisfy P g (k) ^ b 2 P,„(k) and 
a similar relation can be written in terms of b h (which depends on 
halo mass and time) for the dark-matter haloes. 

The introduction of primordial non-Gaussianity has multiple 
observable consequences for the two-point statistics. First, a small 
modification of the matter power spectrum appears on small scales 
(Taruya et al. 2008; Pillepich et al. 2010), due to corrections coming 
from the linear matter bispectrum, which is non-vanishing in this 
case. Secondly, the biasing law between dark-matter halos and the 
underlying mass-density field is altered, becoming strongly scale- 
dependent in the local and orthogonal cases JDalal et al. |2008| 
20081 ISlosar et al.||2008| [Afshordi & Tolley 



Matarrese & Verde 



2008; Valageas 2010, Giannantonio & Porciani 2010, Schmidt & 
Kamionkowski 2010, Desjacques et al.||2011|l. This is due to the 



coupling between long- and short-wavelength modes of the per- 
turbations, affecting the halo power spectrum on large scales. In 
addition to this effect, there is also a smaller scale-independent 
modification to the bias coming from the non-Gaussian form of the 
mass functions ( Slosar et al. 2008 , Giannantonio & Porciani 2010| >, 
affecting the power spectrum on all scales. The above mentioned 
numerous independent calculations and comparisons with N-body 
simulations have shown that the bias at fixed halo mass can be writ- 
ten as 



where <5£>(/nl) and Ai>(&, Tnl) denote the scale-independent and 
the scale-dependent corrections, respectively. If we introduce a 
weighted variance smoothed by a top-hat filter F R {k) 



^ = f -^k"P (k)F 2 R (k), 



(8) 



then the deviation at any given redshift is for each configuration 
approximately given by (Schmidt & Kamionkowski 2010, Des- 
jacques et al. 201 1 1 

2/nl $c b L 



A£ 1oc (£,/nl) - 
W(k,f NL ) ^ 



a(k) 

a(k) o- 2 R0 

6/nl£ 2 °j-2 

a(k) cr 2 



0-R.-2CT' 



1 



,(9) 



where the prime denotes derivatives w.r.t. R. Here the (linear) La- 
grangian halo bias b L = i(/ NL = 0) + 6b(f NL ) - 1 includes the 
(generally small) scale-independent PNG correction. In the spirit 
of the halo model, we will assume that the galaxy linear bias coef- 
ficient b g is given, on large scales, by a weighted average of the halo 
bias at different masses. Note that, on scales larger than ~ IMpc, 
which we only consider in our analysis, the scale dependence is 
translated unaffected from the halo to the galaxy bias. In detail, 
for any choice of galaxy number density n(M) and halo occupation 
distribution (HoD) of mean number Af ? (M), we can write: 



b g (k)oc b h (k, M) n(M) NAM) dM , 



(10) 



which means that, independently from the HoD model, the be- 
haviour of the galaxy bias is the same as for the halo bias: 



b g (k,fy L ) = Wnl = 0) + 6b g (f NL ) + Ab g (k, /nl). 



(11) 



K B {k, / NL ) = fo(/ NL = 0) + 6b(f NL ) + Ab(k, /nl) . 



(7) 



Therefore, when dealing with galaxy clustering, we will 
rewrite the Lagrangian bias as b L = b - 1 with b a nuisance pa- 
rameter to be marginalised over (the index g is understood). From 
the expressions above, we can see that the asymptotical scale de- 
pendence is oc k~ 2 ,k~ [ , k° for the three configurations respectively. 

Fudge factors 

Finally, a note of caution. The Eqs. |9]( have been derived using the 
peak-background split technique to the halo mass function obtained 
through a Press-Schechter ansatz: i.e. assuming that linear density 
perturbations on a given mass scale collapse into dark-matter haloes 
if d > 6 C . To first approximation, it is often assumed that S c - 1.686 
(independent of halo mass) as expected from the spherical collapse 
model in an Einstein-de Sitter universe. However, more rigorous 
applications of the excursion-set approach show that the correction 
to the halo mass function for primordial non-Gaussianity is more 
complicated than this. Using models of triaxial collapse, |Lam &| 
|Sheth| p009| > found that the critical threshold is mass-dependent 
and tends to V0.7(5 C for large halo masses (where the factor 0.7 
is determined by fitting numerical simulations). Alternatively, ac- 
counting for the fluctuations of the collapse threshold measured in 
N-body simulations by Robertson et al. (2009), gives an effective 
threshold 0.89<5 C ( |Maggiore & PJotto|2010b]ay 

Inspired by this theoretical work, in numerical applications, 
the standard collapse threshold 6 C is sometimes rescaled by a fudge 
factor q; this is a coefficient of order unity which is introduced 
heuristically to improve the agreement of the analytical models for 
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the halo mass function and bias with N-body simulations I Carbone 



eTaT1[2008l |Grossi et al.||2009"l |Pillepich et al.|[2010) [Desja cques 
& Seljak||2010| |Wa gner & Verde 2011). However, there are sev- 
eral ambiguities concerning this factor. First, there is no consen- 
sus on the actual role of q. Some authors use it in the ratio (Eq. 
|A10| in the Appendix [A| between the halo mass functions derived 
from Gaussian and non-Gaussian initial conditions (Pillepich et al. 
|20 1 0| [Giannantonio & Porciani 2010). Using the peak-background 
split to derive the halo bias, this gives Eq. (9) with the replacement 
S c — > q6 c ( Giannantonio & Porciani 2010 1. On the other hand, other 
authors use the correction yfqS c in the halo mass function and qS c 
in the halo bias (Grossi et al. 2009). Furthermore, it is sometimes 
assumed that there are two independent fudge factors for the mass 
functions and the bias (Wagner & Verde 201 1). Secondly, even as- 
suming that q is not mass and redshift dependent, significantly dif- 
ferent values of q (with nearly 30 percent deviations) are required to 
reproduce the N-body results when dark-matter haloes in the simu- 
lations are identified using either the friends-of-friends (FoF) or the 
spherical overdensity (SO) method: this makes unclear which value 
should be actually used for galaxies (Desjacques & Seljak 2010| >. 
Finally, the improved model by Desjacques et al. {201 1 1 that we 
adopt, shows good agreement with N-body simulations without the 
need for any fudge factor, at least for a specific choice of the halo 
finder. Because of these reasons, and because the theoretical moti- 
vation favouring a specific value of q is still lacking, in this analysis 
we will set q = 1 throughout. This choice will affect our models for 
the non-linear mass power spectrum (through the halo mass func- 
tion in the halo model) and for the galaxy power spectrum (through 
both the halo mass function and the scale dependence of the galaxy 
bias coefficient). The impact of q on the weak-lensing observables 
is subtle and somewhat degenerate with other cosmological param- 



eters. In fact, in Eq. l[AT0]l /nl is multiplied by a polynomial in q 
whose coefficients depend on the mass variance. On the other hand, 
galaxy clustering studies will constrain the amplitude of PNG al- 
most entirely from the scale-dependent bias on the largest scales. 
This implies that observations will basically constrain the product 
qfytL an d, at least in this case, it will be trivial to correct our results 
whenever future progress on the fudge factor q will be made. An- 
other possibility would be to use q as a nuisance parameter to be 
marginalised over after assuming a theoretically based prior. Even 
though we will not consider this option in this paper, it is clear that 
this marginalisation would degrade the constraints on / NL by a sim- 
ilar amount to the uncertainty in q. 



3 THE NON-LINEAR REGIME 

In order to take advantage of the extended range of measurements 
which will be available from future surveys, it is desirable to ex- 
tend the analyses as far as possible into the non-linear regime. For 
the galaxy clustering, the presence of non-linear biasing compli- 
cates the issue, so to be conservative, it is reasonable to expect an 
accurate theoretical modelling only up to scales k mm = 0.15/i/Mpc 
at z = 0, even though new theoretical methods are being proposed 
which may extend the accuracy of the predictions to smaller scales, 
such as e.g. by |Simpson et al-H201 1) . However, in the case of weak 
lensing, the calculations are made simpler by the absence of the 
bias, and a reasonably accurate modelling is possible up to the 
smallest scales, if the effect of baryons can be neglected or accu- 
rately modelled (see e.g. |van Daalen et"aL] ( |201 1^ for a quantitative 
analysis of this issue). As in most forecast papers appeared so far 
we will not consider this effect in detail here. 



3.1 Halo model 

We compute the non-linear mass power spectrum using the halo 
model develope d by|Ma & Fry|(|2000};rSeljakU2000"|> based on an 
original idea by |McClelland & Silk| < |1977} . See |Cooray & Sheth| 
( 2002 ) for a review. In this approach all the matter in the Universe 
is assumed to be concentrated in a set of discrete and clustered 
dark-matter halos. The non-linear power spectrum can be written 
as a sum of two terms: 

P m (k, z, / NL ) = Pi (k, z, /nl) + Pi(k, z, /nl). (12) 

The two-halo term P 2 dominates on large scales, and represents the 
correlations between mass concentrations lying in different haloes, 
while the one-halo term P\ describes correlations of particles which 
belong to the same halo. These two terms can be calculated as 

r /» f jp(M,z,£) 

n(M,z, /nl) 

Jo 



AOWnl) 

PiihzJw) 



Pm 



dM 



f n(M. 
Jo 

xPo(k,z), 



', z, / NL ) KM, z, k, /nl) P(M,Z,k) dM 

Pm 



(13) 



and there are three ingredients which are needed to implement this 
model: the halo mass function n, the linear halo bias b, and the halo 
density profile p of mean value p,„. The tree-level matter power 
spectrum P is defined in Eq. (|6j; it depends on cosmology through 
its slope, the growth function, and the transfer function, and it does 
not include any /nl correction. There are many possible ways to 
parameterise these ingredients and we have explored different al- 
ternatives trying to maximise the agreement with the power spec- 
tra measured from the N-body simulations presented in Pillepich 
et al. (2010) (hereafter PPH08). Concerning the Gaussian part of 
the halo mass function, we have tried the fitting formula by Tinker 
|et al.| J2008 ), which is based on dark-matter haloes identified with 
a SO halo finder in N-body simulations, and the similar function 
by PPH08, which makes use of a FoF halo finder. For the value of 
the mean overdensity A v which is enclosed by the virial radius of 
haloes, we tried the redshift scaling relation by Brya n & Norman| 
(1998) and a constant value of Av = 200. We also explored the 
halo exclusion prescriptions (Zheng et al.||2002| |Magliocchetti &| 
|Porciani|2003"{ |Tinker et al.|2005) . We found that the best agree- 
ment with the simulations was recovered using the mass function 
by PPH08 and A v = 200, which we will use in the analysis. The 
agreement is at the ~ 10% level as can be seen in Fig. |T| For non- 
Gaussian models, we will then apply the appropriate corrections 
to the PPH08 halo mass function using the method developed by 
|LoVerde^T al. ( 2008 1 and summarised in Eq. jATO}. We will consis- 
tently derive the non-Gaussian halo bias using the peak-background 
split formalism, as reviewed in the Appendix |A| 

We have assumed the analytical description by Navarro et al. 
( |1996| |1997| > for the halo density profile of dark matter haloes 
(NFW). In this approach, the radial density p is written in real space 



(14) 



and it is fully specified by the parameters r s ,p s , describing a scal- 
ing radius where the density profile changes its slope and its associ- 
ated density. Alternatively, we can use the concentration c = Rv/r s , 
where the virial radius R v is defined as the radius of the sphere 
whose mean density is A v = 200 times the average density of the 
Universe. We can finally write the density profile in Fourier space 
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as i |Scoccimarro et al.|200l[|Rudd et al.|2008l > 

p(M, z, k) = 4np s r 3 s I sin(Jfcr,) [Si(/fcr,(l + c)) - Si(fe-,)] 



80, z = 



(15) 



sin(ckr s ) 
kr s (l + c) 



+ cos(/tr,) [Ci(kr s (l + c)) - Ci(fo- S )] 



Here Si(x), Ci(x) are the sine and cosine integrals respectively. The 
concentration is known to depend on halo mass and redshift. We 
use the formula by K lypin et al.| |201 1 ), obtained from the Bolshoi 
simulation: 



c(M,z) = 9.2 k(z)D 1 \z) 
x 1+0.013 



M 



10 12 /j-'M e 
M 



lO^/r'M, 



(16) 



where the function k{z) is introduced to correct for the different 
definition of the virial radius, and it has values k(z = 0) = 1.26, 
and k(z > 1) - 0.96. We have also checked that using the result by 
Bullock et al. (2001 ) for the concentration does not change signif- 
icantly the results. It is worth noticing that the concentration is in 
principle dependent on the PNG parameter / NL but, as was shown 
e.g. by |Smith et aE1p0TT| l; [l)'Aloisio & Natarajan|p0TT) , the ef- 
fect is very small (at least for local PNG) and we will ignore it in 
the rest of this analysis. 



3.2 Power spectrum and comparison with simulations 

With the ingredients described above, we can now calculate the 
full non-linear matter power spectrum P m (k, z); to recover the linear 
regime in the large-scale limit, it is fundamental that 



f n(M, z, /nl) KM, z, k, / NL ) — dM : 

Jo Pm 



1 



(k -> 0) . 



(17) 



We will enforce this constraint in our calculations by adding a con- 
stant to the bias in the smallest mass bin, as also described by Fedcli 
|& Mos cardini (20101. This condition should be automatically sat- 
isfied, but the need to enforce it explicitly arises due to the binning, 
the finite integration limits which exist in practice, and other nu- 
merical issues. 

We then compare the results of the non-linear power spec- 
trum obtained with these recipes with the A-body simulations by 
PPH08. Consistently with PPH08, for our calculations we assume 
the WMAP5 flat ACDM model, with parameters h = 0.701, 
o-g = 0.817, n s = 0.96, Q„, = 0.279, Cl b = 0.0462, £1 A = 0.721. We 
show the result of the comparisons in Fig.[T]at fixed redshift z = 0, 
as the comparison results are very similar at different redshifts. In 
the top panel we can see the full matter power spectra for the choice 
of local /nl = 80, compared with the simulations. Below we show 
the residuals P m (k)/P slmullltlons (k). We can see that the agreement is 
at the 10% level up to scales of 10 h/Mpc. This is the level of ac- 
curacy to be expected from the halo model (Cooray & Shet h|2002| >. 
We can also see a comparison with the method to calculate the non- 
linear power spectrum using a fit from simulations by Smith et al. 
p007] l, whose accuracy is of similar order for k < 2/7/Mpc, but 
then degrades further to the ~ 20% level. 

In Section |6lj| we will compare the results of the Fisher matrix 
analysis obtained using different models for the non-linear evolu- 
tion of the matter power spectrum, finding that the level of uncer- 
tainty still present in the halo model does not significantly impact 
the results of our forecasts in the galaxy clustering case. 
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Figure 1. Comparison of the theoretical matter power spectra with the mea- 
surements from A'-body simulations by PPH08, with local /nl = 80. The 
top panel shows the total P m (k) for the linear (red) and two non-linear mod- 
els: the halo model (blue) and the model by Smith et al. ('2007 1 (green). 
Below we plot the ratio P m (/c)/P s i mu i a tions(^)- All curves are at z = 0: a 
similar (~ 10%) level of accuracy is obtained at other redshifts. 



4 FUTURE GALAXY SURVEYS 

Wide galaxy surveys stand amongst the most powerful ways of 
probing cosmology; they can be used to measure several observ- 
ables, such as tracers of the large-scale structure of the Universe 
and indicators of the expansion rate. For example the Sloan Digi- 
tal Sky Survey (SDSS) ( |Aihara et al.ptJTT] ) and its follow -ups have 
greatly improved our understanding of the standard model, espe- 
cially thanks to the measurement of galaxy clustering and baryon 
acoustic oscillation s (see e.g. |Percival et al.||2007| >, but also us- 
ing galaxy clusters (Koester et al. 20071 and Type la Supernovae 
( |Kessler et al.|2009||L .ampeitl et al. 2009 1. Weak lensing surveys 
(see e.g. |Fu et al.||2008| >" are becoming increasingly powerful in 
constraining cosmology due to the improved control of systemat- 
ics. The same surveys also provide additional cosmological probes, 
e.g. their external correlations with the CMB anisotropies l |Gian-| 
|nantonio et al.|2008| >. 

All of these observables share part of the information, and will 
therefore not be fully independent. A considerable challenge for 
future surveys will be to properly account for their covariance, as 
was discussed e.g. for lensing and cluste ring by|Hu & Jain ( 2004) 
and for lensing and galaxy clusters by |Takada & Bridle| (2007). 
In this paper we study how the constraining power on cosmology 
will improve in the near future using two examples of upcoming 



surveys, focussing on the Dark Energy Survey (DES) (Abbott et al. 



|2005[> for the DETF stage III and an Euclid-like survey iLaureijs 



|et al.|20 09 ) for stage IV. It is reasonable to expect that other surveys 
in the same class will have comparable performances, such as in 
particular the American proposed WFIRST class IV mission. 
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4.1 Stage III: the Dark Energy Survey 

The Dark Energy Survey (DES J] {Abbott et aT1|2005f is an op- 
tical and near-infrared survey which is currently being deployed 
at the Cerro Tololo observatory in the Chilean Andes, and has re- 
cently come to light in late 201 1. The wide-field survey will cover 
5,000 sq. deg. in the Southern sky, reaching ~ 24 magnitude in 
the SDSS bands g, r, i, z and Y. Additionally, the J, H, K bands for 
the same fields are expected to be added from the ESO Vista sur- 
vey. About 300 million galaxies are expected to be observed with 
shapes, photo-z's and positions, as well as 100,000 galaxy clusters 
and 1,000 Type la SNe using a smaller repeated imaging survey. 

In this paper we will study the Fisher-matrix forecasted errors 
on dark energy and PNG using the weak lensing, galaxy clustering, 
and the combined data from this survey. For this purpose, we will 
assume the specifications summarised in Table [T] and that the red- 
shift distribution of the sources can be well approximated by the 
law (Small etal.|1994| > 



dN, 
~dz 



(z) : 



exp 



(18) 



where the parameters are set to a = 2, p = 1.5, in which case 
z is related to the median redshift by zo — z/1.41. We will fi- 
nally split this distribution in eight equally populated redshift bins. 
Due to the uncertainty in the photo-z determination, we will write 
the theoretical redshift distribution of the sources as the convolu- 
tion of these redshift bins with a Gaussian characterised by a dis- 
persion matching the photometric redshift uncertainty cr z (z). This 
approach assumes that the photometric redshift estimation is per- 
fectly calibrated to a spectroscopic sample; the distribution of the 
redshift errors is Gaussian and the r.m.s. error o%(z) is known; catas- 
trophic errors can be identified and downweighted in the analysis. 
See |Kitching etHl] (2008) ; |Bernstein & Huterer| ( |2010fr for more 
detailed approaches which take into account more of the possible 
systematics arising in this case. 

Finally, in the absence of a detailed model we will assume the 
fiducial value of the galaxy bias to follow the law b(z) = Vl + z 
following Rassat et al. ( 2008 1, but we will study the effect of mod- 
ifying the fiducial bias in Section|6!4| 



4.2 Stage IV: Euclid 

Euclid^ is a future space mission of the European Space Agency 
(ESA), which is expected to survey the whole extragalactic sky 
(up to 20,000 sq. deg.) from the L2 point in space (Laurei js" et al.| 
2009). Launch is currently scheduled in 2018. In our analysis we 
approximate this planned survey with the settings and specifica- 
tions described in the Euclid Assessment Study Report (Laureijs 
|et al.|[2009). Although the s pecifications of the survey have since 
evolved ( jLaureijs et al.|201l) , this does not affect the current study 
because our intention is to draw broad conclusions that will be rel- 
evant to the whole class of stage IV surveys including WFIRST, 
LSST, and others. See in any case Section [6.10| for a discussion of 
how the results change when using the latest Red Book specifica- 
tions. The Euclid mission is expected to perform two main surveys, 
photometric and spectroscopic. 

The photometric part should measure photo-z and elliptici- 
ties in the optical and near-infrared bands (one broad visual band 
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Figure 2. Redshift distributions used for the forecasts of both photometric 
and spectroscopic data sets for the Euclid satellite and the DES. The photo-z 
distributions are given by an analytic function jSmail et al.|1994) , while the 
spectroscopic part is numerically estimated by Geach et al. (2010), where 
the Euclid specified flux cut is used, 4 ■ 10~' 6 erg s _1 cm -2 . The distribu- 
tions have been already convolved with the probability density function of 
redshift measurement errors. 



R+I+Z and Y, J, H IR bands), up to mag ~ 24.5 in the visual and 
24 in the IR. The requirement specifications are described in Table 
[T] The expected number of observed galaxies is of the order of a 
billion. For this survey we will also use the approximated redshift 
distribution by |Smail et al.|ft994) , dividing the sample in 12 red- 
shift bins, whose distribution is again convolved with the expected 
photometric errors, as shown in Fig. [2] The photometric galaxies 
are distributed in the redshift range < z < 2.5. 

The spectroscopic survey is expected to use a slitless spec- 
trometer which will mainly detect the Ha emission line of galaxies. 
The spectrometer will have a resolution A/AA = 500, giving a red- 
shift uncertainty of cr z (z) = 0.001(1 + z). The wavelength range of 
this instrument will be limited to 1000 nm < A < 2000 nm, meaning 
that only galaxies at 0.5 < z < 2 will have measurable Ha lines and 
thus redshifts. The limiting flux is placed at 4 • 10~ 16 erg s -1 cirT 2 , 
which combined with the expected success rate of the spectrometer 
e = 35% yields ~ 60 million galaxies, using the predicted tabulated 
calculations by Geach et al. ( 2010), which was based on empirical 
data of the luminosity function of Ho- emitter galaxies out to z = 2. 
We will use this tabulated prediction as our fiducial redshift dis- 
tribution of the sources, consistently with the Euclid Study Report 
specifications (Laureijs et al. 2009). The remaining specifications 
are again to be found in Table [T] We will finally split the distribu- 
tion into 12 equally populated redshift bins, as we can see in Fig. [2] 
Again, we will take the fiducial val ue of the galactic b ias to fol- 



Rassat et al. 



(2008), which is 



http : / /www . darkenergy survey . org/ 
http://sci.esa.int/euclicl 



low the law b(z) = Vl + z following ' 
a good approximation to recent studies from semianalytic models 
of galaxy formation such as e.g. by jOrsi et al.| (2010), but we will 



study the effect of changing this choice in Section|tx4| 
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Parameter 


description 


Euclid photometric 


Euclid spectroscopic 


DES 


<Tz(Z)l\i + Z) 


redshift uncertainty 


U.UJ 


n nni 


u. i 


z 


median redshitt 


1 


j ] 


U.O 


n 


galaxy density 


40 arcmin -2 


0.7775 arcmin 


1 2 arcmin 


A 


surveyed area 


20 000 sq deg 


20 000 sq deg 


5 000 sq deg 


y 


intrinsic ellipticity noise 


0.247 




0.16 


dN/dz(z) 


galaxy distribution 


Smail et al. 


Geach et al. 


Smail et al. 




parameters of gal. dist. 


a = 2 


flux cut = 4 ■ 10~ 16 erg/s/cm 2 


a = 2 






(3 = 1.5 


efficiency = 35% 


fi = 1.5 


M 


number of redshift bins 


12 


12 


8 



Table 1. Technical specifications of the surveys used for forecasts. 



5 FISHER MATRIX FORECASTS 



5.1 Formalism 



Let us assume, following e.g. Teg mark et al.| (19971; Tegmark 
{T997jl, that we have a data vector x, whose elements are random 
variables which depend on the cosmological parameters 1?, with a 
likelihood function L(x, Then, if we fix X = - InL, the Fisher 
matrix can be defined as the ensemble average 



/ d 2 £ 



(19) 



and the Cramer-Rao inequality proves that its inverse represents the 
best possible covariance matrix for the measurement of the param- 
eters 1? with the help of unbiased estimators. In the limit of large 
data sets, if the distribution becomes Gaussian, this inequality be- 
comes an equality, and the standard deviation on a given parameter 
??, is given by <x(j?,) = V(F~ '),■,-. Under the assumption of a Gaus- 
sian likelihood function L, for data of mean fj, = (x) and covariance 
matrix £ = <xx r > - fifi T 
sum of two pieces: 



the Fisher matrix can be written as the 



F w = iTr[s- 1 S,, S- 1 S, J ] + At J £" 



(20) 



where the commas denote derivatives with respect to the parame- 
ters. 

To calculate the Fisher matrix, we need first to define our pa- 
rameter set and its fiducial model, which we take in our base anal- 
ysis to be for the galaxy clustering case: i9 = {^a = 0.721, £l b = 
0.0462, Sl m = 0.279, h = 0.701, n„ = 0.96, cr 8 = 0.817, w = 
-1.0, w a = 0, /nl = 0, b = bfid), where the usual set of cosmologi- 
cal parameters is chosen to match the WMAP5 set used in Pillepich 
|et al.|p010| , the dark energy equation of state is parameterised as 
w(z) = w + w a z/(l + z) (Chevallier & Polarski 2001 1, and the array 
b of fiducial value b fld is a set of nuisance parameters to describe 
the total scale-independent part of the galaxy bias in each redshift 
bin (see discussion in Section |2T2) . Notice that the degeneracy be- 
tween o" 8 and b is broken by the non-linear corrections to the power 
spectrum at k > O.IA/Mpc. For the lensing case, as there is no bi- 
asing, we do not include bias nuisance parameters. We will also 
consider different parameter arrays and extra nuisance parameters 
in some cases, as described in the following sections. 



5.2 Observables 

Starting from a 3D random field (p(x), its 3D power spectrum P{k) 
can be defined as in Eq. {2}. Similarly, if we consider a 2D ran- 
dom field on a sphere £ (n), the angular power spectrum C\ is given 
by {aim at ,) = 6u>S mm > Ci, where a tm are the spherical harmonic 



coefficients of the field. We will use in our analysis the informa- 
tion contained in all the relevant two-point correlation functions 
(or power spectra) which are measurable by the the DES and Eu- 
clid missions: we shall consider the 2D power spectra of weak lens- 
ing (using the photometric redshift catalogue) and galaxy clustering 
(using both the photometric and spectroscopic catalogues), and the 
3D power spectrum for clustering (spectroscopic), as described be- 
low. We will first present the forecasts for each single probe, and 
finally combine lensing and clustering, including all the relevant 
two-point correlations and covariances. For all numerical calcula- 
tions we use the iCosmo package by Refregier et al. ( 201 1) . 

An important point is whether this basic Fisher setting, which 
assumes a Gaussian distribution of the covariance, is accurate 
enough. More extended calculations have shown that extra terms 
arise, due to the connected matter trispectrum and also to the cou- 
pling of large-scale modes which are outside the volume of the sur- 
vey ( |Sato et al.|[2009| |Takada & JaTn]|2009"l >. The effect of these 
terms is to reduce the signal-to-noise ratio obtained from the small 
scales, thus mostly affecting the cosmic shear results. The worsen- 
ing of the lensing marginalised errors is expected in most cases not 
to exceed the ~ 10% level (Tak ada & Jain|20 09 1 and thus, while it 
will be important for data analysis, we will ignore it at the forecast- 
ing stage. 



5.3 Cosmic shear 

We can describe a tomographic galaxy survey of iV galaxies in M 
redshift bins with a set of distribution functions for the galaxy den- 
sity fluctuations in each f-th bin dNj/dz(z), so that the total number 
of galaxies in the bin is 



r 00 , dN, , 

N '= *T7fe>< 
Jo az' 



N = Y J N t . 



(21) 



In a Friedmann-Robertson- Walker universe of total energy density 
Q.0 we can write the curvature of the hypersurfaces at constant cos- 
mic time as 



Hi 

K= -£(£2o-l), 



and the comoving angular diameter distance is 

( K- [ > 2 sin(^ 1/2 r) if£>0 

r K (r) = \ r if K = , 

( (-K)- 1 ' 2 sinh(-K 1/2 r) if AT < 



(22) 



(23) 



which coincides with the comoving distance r in the flat case. The 
spatial part of scalar metric perturbations can be decomposed in 
terms of the solutions of the Helmholtz differential equation {Ab-| 
|bott&Schaefer|1986| > 
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(D 2 +k 2 ) <2(r) = 0, 



(24) 



where D 2 is the covariant Laplacian in curved space. The eigen- 
functions Q(r) generalise the concept of plane waves to curved 
spacetime. Introducing spherical harmonics leads to variable sep- 
aration and the radial part of the eigenfu nctions can be written as 



R kl (r) = O^(r) (Abbott & Schaefer 



1986), y?= ik 2 + K, and *f (r) 
are the ultra-spherical Bessel functions (note that Abb ott & Schae-| 
|fer| (T986 ) write their explicit form in terms of a different radial 
coordinate). 

We can use the weighting kernels for cosmic shear ( |B artel- 1 
|mann & Schneider|2001) 



W € >(z) 



1 



3H Q,„ 



N, 2c 2 a(z) 



, ,dN r K [r(z')-r(z)] 
dz' r K [r(z')] 



(25) 



to calculate the projected 2D cosmic shear power spectrum between 
the bins (, j as 



C 



- f dpp 2 f dnW(n)tf(n) 

n Jo Jo 

• f rfr 2 W f '(r 2 )*?(r 2 ) • P m (J3,r u r 2 ), 
Jo 



(26) 



where P „, here denotes the cross-spectrum between the matter dis- 
tribution at two different cosmic epochs corresponding, on our past 
light cone, to distances r t and r 2 .[^]This can then be written in the 
Limber approximation as 



f 



dr 



W e '[z(r)]W e i[z(r)]P r , 



1 + 



r K (r) 



;z(r) 



(27) 



and the full 3D matter power spectrum can be calculated e.g. using 
the halo model formalism of Section[3] 

The observed lensing power spectra C ( ' 1 may be modelled 
as a sum of the theoretical spectra and the noise. For each pair of 
redshift bins i, j, we have 



Cp=Cp+Np, 

where the noise term N? ej is due to the shape noise 



AC 



(28) 



(29) 



Here n, is the angular number density in the i-th bin, and y is the 
r.m.s. shear arising from the intrinsic ellipticities of the galaxies, 
which is a specification of the survey (see Table[TJ. 

Another important source or errors in cosmic shear are intrin- 
sic alignments: we do expect a fraction of neighbouring galaxies 
will have correlated ellipticities due to tidal fields. This effect in- 
troduces a systematic into the cosmic shear analysis which has to 
be taken into account when constraining cosmology from real data 
( |Joachimi & Bridle|2010||Kirk et al.|20T0) , while we will neglect 
it for our forecast as the biasing introduced on the parameter esti- 
mation typically does not exceed the ~ 10% level. 

We can now calculate the Fisher matrix element from Eq. \2Q) 
for the case of weak lensing. A common assumption in this case is 
to identify the observables x with the afj instead of the C^' , be- 
cause the a h „ are at least in the linear regime Gaussian. This brings 
in the advantage of cancelling the second term of Eq. l |20| l, since 
Mi = ( a \ m ) = 0, Vj;, /, m). By definition, the variance of each of 



the (21 + 1) is given by YL ab = C^'Sai,, and thus Eq. |20| yields 
l |Hu & Jain|2Ti04"l|Arnara & Refregier|2007l > 



: /sky X 

'='mir 



(21 + 1) , 



(30) 



where / sky is the sky coverage of the survey. Here C ; E is a matrix of 
dimensions M x M, whose elements are the lensing power spectra 
at a fixed / between each pair of redshift bins C ; f|£ ' . Finally, D e ff is a 
matrix containing the derivatives of the spectra with respect to each 
cosmological parameter ft a , whose elements are defined as 



D, 



dC, 



(31) 



We shall calculate the forecasted weak lensing power spectra ac- 
cordingly with the proposed specifications of the DES and Euclid 
missions, as described above. 

The range of multipoles used in the sum of Eq. < |30| > [/ m in, /max] 
is a sensitive issue, since we would like to set it as broad as possible 
to use all available information. The value of / max is related to how 
far we can trust our modelling of the non-linear regime. In the case 
of weak lensing, as we are free from the problems of non-linear 
biasing, we take / max = 20,000, which gives = 8.5/;/Mpc at 
z = 1 assuming the best-fit WMAP cosmology. Again this assumes 
that the effect of baryons to the matter power spectrum can be either 
ignored or accurately modelled. Another issue which may affect the 
analysis at these small scales is the non-Gaussian contribution to 
the covariance; we do not account for these corrections whose ef- 
fects are expected to be small, as discussed at the end of Section |5"72] 
Further, the theory calculation in this regime is strongly dependent 
on the parameters of the halo model, such as the concentration, the 
subhalo distribution etc. We will discuss in Section lo77l the effects 
of changing this limit, where we will find that with more conserva- 
tive choices the differences are still small (see Fig.|7|: for example, 
the marginalised error on / NL including Planck priors degrades in 
this case by only 15% between / max = 20, 000 and / max = 5, 000. 

The minimum multipole used is also important. It is expected 
that the accuracy of the Limber approximation will deteriorate for 
/ — > 0. Here, to be complete and conservative, we will only use this 
approximation at small scales, for / > 200, and will use the com- 
plete exact formula at larger scales (using the linear approximation 
for the cross spectrum at two times). In this way, we can extend the 
calculation all the way to the largest scales, and will take / mill = 5. 
However, as shown in Fig. [3] the exact calculation is important only 
for the case of galaxy clustering, where it can not be neglected as 
in the presence of non-Gaussianity, the large scales (small multi- 
poles) are the most affected by the scale-dependent bias. Note that 
the Limber approximation overestimates the power on the largest 
scales for positive / NL while it underestimates it for a Gaussian 
model. Since the Fisher matrix is computed taking the derivatives 
of the signal with respect to model parameters, the Limber approx- 
imation is rather inaccurate for a fiducial with / NL = 0. 

We can see the resulting 2D marginalised forecasts for a 
Euclid-like survey for local PNG in Fig. [4] and all the marginalised 
ID error bars in Fig. [6] The results are also summarised in detail 
in Tables [2] [3] [4] for the three PNG configurations. We present the 
results for DES in Table|5]for the local case. 



5.4 Galaxy clustering, 3D 

The second observable we c< 

> /. which is typically measured via the 3D power spectrum using the 



3 For a universe with positive curvature the integral over the corrected 

wavenumber yS should be replaced by a discrete sum such that j} > 3 and The second observable we consider is the clustering of galaxies 
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Figure 3. Summary of all observables used in this work for the fiducial 
model described in Section|5] In the upper panel, we show the projected 2D 
spectra for lensing, galaxy clustering, and their cross-spectra for photomet- 
ric and spectroscopic surveys. In the bottom panel we show the 3D galaxy 
power spectrum. Note that the Limber approximation is inaccurate on large 
scales for the galaxy-galaxy spectra, as well as the linear-theory power 
spectrum is on the small scales (dashed lines). For the lensing and galaxy- 
lensing cases, the Limber approximation works well due to the wider distri- 
bution of the sources. The shaded areas represent cosmic-variance errors for 
a half-sky survey. A redshift bin centred around z = 1 is used in all cases. 



spectroscopic part of the survey. In this case we then consider 
measurements of the galaxy power spectrum for each z-th red- 
shift bin P g .(k) = b 2 g .(k)P m (k), where b g .(k) = \b gt +Abi(kj\ is 
the galaxy bias of the ;-th bin. The scale-dependent part of the 
bias is computed as explained in Section [2] while we consider the 
scale-independent part as a nuisance parameter to be marginalised 
over. Due to the high accuracy of the spectroscopic redshifts, we 
can ignore any cross-spectra between different redshift bins, since 
the redshift distribution functions will have negligible overlap, and 
we will also neglect the effect of covariances due to the longitudi- 
nal modes. We also include corrections due to two effects, as de- 
scribed e.g. by |Seo & Eisenstein|p0fJ3) ; |Song & Percival| < [2009) : 
the redshift-space distortions and the Alcock-Paczynski effect. 

Redshift-space distortions are caused by the peculiar motion 
of galaxies ( Kaiser 1987 1, and produce an enhancement of the over- 
density field 6 proportional to (1 + ffi 2 ), where / is the logarithmic 
derivative of the growth function with respect to the expansion fac- 
tor a, and /j is the cosine of the angle to the line of sight. If we 



can exclude velocity bias, the correction of the power spectrum be 
comes (see e.g. |Peacock & D odds ( 1994}) 



J*(k) = J»„(k) \\+tfp m m\ f 



MHz) 



(32) 



where fi gl {k) = f/b gi (k), cr,., is the ID pairwise velocity dispersion 
of the galaxies, and H(z) is the Hubble function. We will use for 
the function F the form F(x) = e~ x (Pe rcival et al.| [2007). This 
gives a small correction whose cosmological content can be ignored 
and treated as a nuisance parameter to be marginalised over ( |Song| 
|& Percival|2009] >. In more detail, we can decompose the velocity 
dispersion in each redshift bin i as the sum in quadrature of an 
intrinsic term coming from the finger-of-god (FoG) effect, o"f g,/> 
and a term due to the redshift uncertainty of the survey cr. given in 
Table [T| as 

2 

2 /i \2 ^roG,i 2 2 

= (1 +zY + c z ai . (33) 

Since the redshift evolution of this intrinsic term is uncertain, it 
would be desirable to introduce independent parameters <tf g,; for 
each redshift bin, and then marginalise over them. However, we 
found that, given that this factor enters the calculation together with 
the function l/H(z) in Eq. {i2\ , leaving these parameters free to 
change brings in severe degeneracies with any other cosmological 
parameter which is able to alter H(z). For this reason, we have been 
forced to adopt a less conservative approach, where we choose a 
functional form 



+ z 



(34) 



and only marginalise on one parameter <tf g,o- This scaling with 
redshift assumes that the observed galaxies reside in halos with 
nearly constant mass. It is also possible to make this approach more 
conservative, e.g. by considering more complicated parameterisa- 
tions for this function, but again the constraints on all parameters 
which enter the Hubble expansion H(z) will worsen dramatically. 
As a fiducial value, we take ctfog.o = 250 km/s, as the Euclid spec- 
troscopic sample is made of star-forming galaxies, which are not 
generally located in massive haloes, and are thus expected to have 
a low velocity dispersion. 

On the other hand, the Alcock-Paczynski effect occurs since 
to infer galaxy distances from the observed redshifts and positions 
we have to use a reference cosmology, which is different from the 
'true' one. This can be corrected, and the true values are ( |Seo &| 
|Eisenstein|2003| > 



Kl (1 - M ^ —-, + ( fc ref ^ref) 



' Dljz) 



HIM 



true V 

l 



KM 



1/2 



KM K 



(35) 



Here we will identify the reference cosmology with our fiducial 
model for simplicity. Finally, the observable corrected power spec- 
trum, which we will use to calculate the Fisher matrix, is given by 



PgAktmet Mtrue) 



DLMHrAz) 
Dl f (z)H lme (z) 



x P m (k tme ) F 



^(*true)[l +A#trae) 

^trueA'l 



2 

A* true 



(^true£ i true C v,/ 1 
H\z) ) 



+ P, 



shot,g ; ■ 



(36) 



The last term on the r.h.s. is due to shot noise. If this is Poisso- 
nian, its fiducial value is P s hot,a = However, it is possible to 
be more conservative and assume that the mean galaxy density is 
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not perfectly known, or that there are other sources of shot noise. 
We will include these unknown contributions by introducing one 
additional nuisance parameter for it in each bin, and marginalising 
over them. 

The Fisher matrix calculation is here based on the works by 
Feldman et al. ( 1994); Tegmark ( 1997) amongst others. In the 3D 
case, it is customary to identify each of the observables x, with the 
average power in a thin shell of radius kj in Fourier space, of width 
dkj and volume V, = Ank 1 i dk i l(2nf . In this case, for each redshift 
bin a and angle yi b we have non-zero means /i, = P gci (kj,ji b ), and 
covariances 



V,V eff (£,,jUj,) 

where the effective volume is 
h{r)P ga (k,n) 



d 3 r, 



(37) 



(38) 



and n{r) is the selection function of the survey. The dominant term 
of the Fisher matrix is now the second one in Eq. \20) which, adding 
the effects of redshift-space distortions, can be finally written as 
( |Seo & Eisenstein|2003l 



din P(k,fi) din P(k, fx) 



nk 1 



-dkdfi. (39) 



7, mm W a d&f) c " v """' (2nf 

Since here we are not accounting for bias non-linearities nor scale- 
dependences at small scales, our calculated galaxy power spectrum 
is exact only in the linear regime; therefore to limit the effects of 
non-linearities, in particular for what concerns galaxy biasing (see 
|Roth & Porciani (2011)), we cut our calculation at a scale £ max = 
0.15 /i/Mpc at z = 0. We can then make this limit more meaningful 
by evolving it as a function of z, by imposing the condition 

^- i k 2 P (k,z) = const.. (40) 

Here we have set k min = 10~ 3 ft/Mpc, as this roughly corresponds to 
the effective volume of our Euclid-like survey. The obtained values 
of fcmax(z) can then be translated into / m;lx by using the approxima- 
tion 



max max 



(41) 



In this way the choice £ raax = 0.15 /?/Mpc at z = translates into 
imposing that the variance is fixed to tr 2 (z) = 0.36, i.e. the r.m.s. 
of the perturbations is cr m 0.6, meaning that perturbation theory 
is generally respected (peaks of order S ~ 1 are rare), and so this 
is a conservative choice. To be even more conservative, we impose 
that in any case the smallest scale can not go above k^ x = 0.3 
/z/Mpc, as was similarly done by Wang et al. (2010). In practice, 
this maximum k is reached around z — 1. We have checked that 
at this redshift, the non-linear contribution to the power spectrum 
does not exceed ~ 13%, and is decreasing at higher redshifts. 

We can see the results for Euclid, local PNG in Fig. [4] and all 
the marginalised ID error bars in Fig. [6] The results are also sum- 
marised in detail in Tables|2][3]|4]for the three PNG configurations. 
We remind the reader that in this case the results are marginalised 
over nuisance parameters for the Gaussian part of galaxy bias, the 
Finger-of-God effect, and where indicated also shot noise. 

5.5 Galaxy clustering, 2D 

In addition to the 3D galaxy clustering in redshift space discussed 
in the previous section, we have also calculated forecasts for the 



projected clustering on the celestial sphere. This is done for sev- 
eral reasons: first, in this way it is much easier to combine the re- 
sults with weak lensing in a consistent way, as we will discuss in 
the next section. Secondly, this allows us to directly compare the 
performances of the photometric and spectroscopic surveys of the 
Euclid-like mission in terms of galaxy clustering: a priori, it is not 
obvious whether the tighter constraints on cosmological parame- 
ters could come from the more numerous galaxies with photometric 
redshifts or from the less rich but more accurately located spectro- 
scopic sample. Thirdly, we can thus present forecasts for the galaxy 
clustering of the DES, which does not include, at least in its initial 
form, a spectroscopic survey. 

The 2D galaxy spectrum between a pair of redshift bins i, j is 
a projection of the 3D spectrum which, in analogy with the weak 
lensing case above, can be written as 

cf Si = - f dpp 1 f AiWCftrOoffa) 

n Jo Jo 

■ f dr 2 W'(0, n) ®f(r 2 ) • PJfi, n , r 2 ) ■ (42) 
Jo 

This can then be written in Limber approximation in analogy with 
Eq. \21\. The sources in a bin i are now given by 



Nj dr 



(43) 



The observed 2D galaxy power spectra C, 1 ' can again be modelled 
as a sum of the theoretical spectra and the noise. For each pair of 
redshift bins i, j, we have 



■SiSj 



C, J = C,' + N\ 
where the noise term N, ' 1 is now due to the shot noise: 



AC 



1 



(44) 



(45) 



The Fisher matrix in this case is defined exactly as in Eq. l |30[ ) 
above, with the lensing power spectra replaced by the galaxy- 
galaxy ones. 

For completeness and for the purpose of comparison, we will 
use in this case both photometric and spectroscopic data sets for 
Euclid, as described in the above sections. Note that, as already 
mentioned for the 3D case, since here we are not accounting for 
bias non-linearities nor scale-dependences at small scales, Eq. l |42[ > 
is exact only in the linear regime; therefore we use the same limit 
on the small scales described above for both cases, i.e. we assume 
now k max = 0.15 /i/Mpc at z = 0, evolved in redshift as in the 3D 
case described above. We will discuss in Section loTTl the effects of 
changing this limit. Also in this case, we use the Limber approxi- 
mation for / > 200, and the exact calculation for larger scales, so 
that we can extend the analysis all the way to Z m j„ = 5. See Fig. [3] 
for a comparison of the Limber power spectra with the exact cal- 
culations for the range of projected observables we consider. Here 
we can see how the Limber approximation departs from the exact 
calculation at large scales for the galaxy spectra. 

We can see the results in Fig. [4] for both photometric and 
spectroscopic parts of Euclid for the local PNG case, and all the 
marginalised ID error bars in Fig. [6] The results are also sum- 
marised in detail in Tables [2] [3] [4] for the three PNG configura- 
tions. These results are marginalised over nuisance parameters for 
the Gaussian part of galaxy bias. 
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5.6 Combined results 

Finally, we shall combine the constraints from lensing and galaxy 
clustering in the projected case, as it is the only case where it is 
straightforward to do it consistently. Since the two probes are based 
on the same density field, it is crucial to correctly include their full 
covariance; on the other hand, it is also possible to include the 
lensing-clustering cross-correlation as a signal. These two opera- 
tions can be achieved by considering one single Fisher matrix, as 
done by Hu & Jain ( 2004 ), and defined by 



(46) 



In this case, the matrices C*, D* will contain all the combinations 
of 2D spectra: C,' ej , C ( e,s \ Cf Sj . Notice that the noise on the cross- 
spectra is zero, since we are assuming that shape and Poisson noise 
are uncorrelated: C ( e,Sj = Cf 1 , 

We consider in this case the combination of the lensing fore- 
casts using the photometric catalogue with both photometric and 
spectroscopic galaxy clustering. The combined results for Euclid 
are reported for both cases in Fig.[5]for local PNG, and summarised 
in detail in Tables [2] [3] [4] for the three PNG configurations. These 
results are marginalised over nuisance parameters for the Gaussian 
part of galaxy bias. 

All the marginalised ID error bars are shown in Fig.[S]for the 
local PNG case. We also present the results obtained for the DES 
specifications in Table [5] for the local case. We can also see from 
Table § the comparison between naively summing the lensing and 
clustering Fisher matrices versus performing the full analysis. We 
can see here that the forecasted errors for the complete combina- 
tion are smaller than the simplistic sum, in average by a factor of 
two. This points us to conclude that the effect of including all the 
cross-spectra between lensing and clustering as observable signals 
in the analysis is more important than the degrading effect which 
is produced by considering the covariances between the clustering 
and lensing auto-correlations. 



6 DISCUSSION AND ADDITIONAL ANALYSES 

We present here some additional results which extend our basic 
analysis of the previous section, and discuss some important com- 
plementary issues. 



6.1 Additional priors 

We first study the effect of adding as a prior the forecasted con- 
straints for the Planck CMB satellite, whose results will be avail- 
able by the time these surveys are completed. The Fisher matrix 
for Planck is obtained in analogy with the galaxy 2D case, and the 
CMB spec tra are calculated with the CAMB CMB code (Lewis 
etal.|2000). We describe this in more detail in the Appendix |B| | 

The reduced likelihood intervals are also shown in 
Figs. [4] [5] [6] We can see that in this case all the constraints im- 
prove greatly, with the significant exception of / NL . This is because 
the Planck prior we are using is only taking into account the CMB 



4 We acknowledge Jochen Weller for supplying us with the CMB Fisher 
matrix code. 



power spectrum, and not the bispectrum, which is where most in- 
formation on non-Gaussianity is. In other words, our Planck Fisher 
matrix does not depend on /nl at all. 



6.2 Summary of the results 

We summarise in Fig.|6]all the marginalised ID constraints which 
are obtained from the Fisher matrix analysis in the Gaussian and 
non-Gaussian cases, with and without the imposition of the Planck 
priors. Here we can also see a comparison between the fore- 
casted future results for Euclid and the current constraints from 
the WMAP satellite and other probes. It is particularly interesting 
to notice that the errors on the standard cosmological parameters 
remain largely unchanged when /nl is added. We also found by 
looking at the covariant elements in the Fisher matrices that _/nl is 
almost completely uncorrelated from the other parameters, which 
is due to the very particular scale-dependent behaviour which it 
produces on the observables. 



6.3 The different shapes of non-Gaussianity 

We now focus on the differences between the three types of primor- 
dial non-Gaussianity described in Section|2] namely the local, equi- 
lateral and orthogonal bispectrum shapes. In all these cases PNG 
alters both the matter power spectrum through the halo mass func- 
tion (which depends on the skewness of the linear density field), 
and the galaxy bias. 

The different forms for the skewness are described in the Ap- 
pendix [A] while we use the expressions for the corrections to the 
bias in Eqs. {9} d erived by [ Schmidt & Kamionkowski (2010} and 
improved by Desjacq ues et al.| ( |2011| >. These expressions are ap- 
proximated but are rather accurate at small wavenumbers for the 
local and orthogonal cases ( Schmidt & Kamionkowski 2010) where 
the scale-dependent part of the bias scales as k~ 2 and k' 1 , respec- 
tively. On the other hand, in the equilateral case, the bias becomes 
asymptotically constant for small k and the only scale-dependence 
of the bias comes from the linear transfer function appearing in the 

function a. Unfortunately, the numerical factor which appears 

in the bias correction is typically large, so that the scale-dependent 
deviations sourced by T(k) easily outweigh the effect of the mass- 
density skewness in the matter power spectrum at large k. This is 
not ideal, as Eqs. (|9j have been derived with the peak-background 
split technique, which should only hold on large scales. For this rea- 
sons we decide here to be conservative, and to keep the bias fully 
scale-independent in the equilateral case, fixing it at all scales to its 
small-i value. 

The marginalised constraints on the non-Gaussianity parame- 
ter /nl are shown in Table^Jfor the three configurations. Here we 
can see that the constraints in the orthogonal and equilateral cases 
get generally weaker for the clustering probes due to the reduced 
importance of the scale-dependent bias. The constraints from lens- 
ing remain nevertheless quite strong, as the effect of the skewness 
is independent from the bias; the results for the combined probes 
are therefore very promising for all the configurations. 

As a further note of caution, we have nevertheless to stress 
that, differently from the local case, we could not check with N- 
body simulations the predictions of the non-linear regime for the 
orthogonal and equilateral cases. This means that the constraints 
from lensing and its combinations in these two cases are necessar- 
ily less reliable than both the local case and the constraints from the 
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Figure 4. Fisher matrix forecasts for the Euclid-like survey, using weak lensing (photometric survey), 2D galaxy clustering (photometric and spectroscopic 
surveys), and 3D galaxy clustering (spectroscopic only). For lensing, the used multipoles are from / ra j n = 5 to Z max = 20 000, while for clustering the maximum 
mode is ir max = 0.15 /i/Mpc at z = 0. The forecasted posteriors are marginalised over the other not shown parameters. The blue ellipses refer to Euclid only, 
while in red we show the results including Planck priors. Notice that the axes ranges are different, which is necessary given the different constraining power 
of the different observables. 



Euclid data 


o-(fi A ) 


o-(fift) 


o-(n m ) 


o-(h) 


o-(n) 


o"(o"s) 


o-(wo) 


o-(w a ) 


o"(/nl) 


lensing 

2D clust. phot. 
2D clust. spec. 
3D clust. 
l-shot noise 


.035 (.014) 
.064 (.037) 
.087 (.019) 
.063 (.0079) 
.076 (.0082) 


.012 (.0016) 
.0085 (.0043) 
.0067 (.0027) 
.0041 (.00096) 
.0044 (.00098) 


.011 (.0081) 
.050 (.026) 
.023 (.015) 
.018 (.0036) 
.022 (.0036) 


.10 (.011) 
.12 (.033) 
.12 (.020) 
.020 (.006) 
.026 (.0062) 


.033 (.0040) 
.027 (.0036) 
.059 (.0039) 
.024 (.0033) 
.032 (.0035) 


.017 (.011) 
.042 (.028) 
.042 (.020) 
.030 (.016) 
.044 (.019) 


.12 (.087) 
.51 (.29) 
.12 (.092) 
.10 (.030) 
.13 (.031) 


.55 (.34) 
1.6 (.94) 
.91 (.33) 
.53 (.090) 
.65 (.094) 


37 (19) 
3.2 (3.0) 
6.7 (6.2) 
4.4 (4.2) 
4.4 (4.2) 


lens. + 2D phot, 
lens. + 2D spec. 


.014 (.0069) 
.017 (.0083) 


.0035 (.00077) 
.0049 (.00094) 


.0034 (.0032) 
.0048 (.0044) 


.033 (.0050) 
.044 (.0064) 


.01 1 (.0030) 
.014 (.0035) 


.0035 (.0032) 
.0051 (.0047) 


.037 (.035) 
.052 (.047) 


.19 (.13) 
.24 (.17) 


2.8 (2.6) 

4.9 (4.6) 



Table 2. Expected marginalised errors from Euclid-like data for all parameters for the local non-Gaussianity case. The numbers within parentheses include 
the forecasted priors from Planck CMB temperature power spectrum. The additional line for the 3D clustering includes marginalisation over the shot noise 
parameters in each redshift bin. 
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Figure 5. Full combination of Euclid-like data. Fisher matrix forecasts for the combination of weak lensing (photometric survey) plus 2D galaxy clustering 
(photometric and spectroscopic surveys). Both the range of used multipoles and the colour coding of the ellipses are as described above for the previous plots. 
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Table 3. Expected marginalised errors from Euclid-like data for all parameters, for the orthogonal configuration. The constraints on /ml from the scale- 
dependent bias are degraded, while the results from lensing are in this case stronger. 



scale-dependent bias, whose predictions are well tested and under- 
stood. 

Finally, we can compare Tables [2] [3] [4] to see how much the 
errors on the other parameters change. Here we can see that when 
the constraint on /nl degrades greatly, such as in the galaxy clus- 
tering probes alone, the constraints on some other parameters also 
degrade significantly, especially erg. However the combined lens- 
ing+clustering result stays largely unchanged. 

Our constraints on / NL are weaker than what was found by 
|Fedelietal.|pOTT) for the 3D case only, as these authors calculated 
forecasts by varying only erg and /nl while fixing the rest of the 
parameters to ACDM and also used a more optimistic cutoff than 
ours for the non-linear regime. 

6.4 Dependence on the fiducial 

We then study the dependence of our forecasts on the choice of 
the fiducial values for the parameters, in particular extending the 
analysis to / NL J= (in the case of local PNG only), or w t -1. For 
this purpose, we run our full analysis for two non-standard cases, 
where the fiducials have the value of one parameter altered to wo = 
-0.95 and / NL = 30 respectively. 

Fiducial with dark energy 

In the case of dark energy, we calculated the forecasted errors 
changing the value of the fiducial equation of state from w = - 1 
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2.8 (2.6) 


2.9(1.6) 


12 (6.5) 
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4.9 (4.6) 


3.4(1.9) 


11 (7.0) 



Table 7. Forecasted results for non-Gaussianity for the different configu- 
rations. Constraints from galaxy clustering weaken in the orthogonal and 
equilateral cases, as the scale dependence of the bias is reduced. For the 
lensing (and thus also combined) cases, the constraints remain strong, as 
the non-Gaussian mass function is in all configurations expected to deviate 
at similar levels from the Gaussian model. 



to w = -0.95. We obtain that the constraints do not change signif- 
icantly, as we can see in Table [8] we can conclude that the model 
with dark energy of wq = -0.95 would be effectively detected by 
Euclid, even without the need of Planck CMB priors. 



Non-Gaussian fiducial 

We chose to look at a local non-Gaussian model with / NL = 30, 
as this is near the peak of the current posterior probability distribu- 
tion given current data. We found that setting the fiducial to such a 
model degrades to some extent the accuracy of the errors. For the 
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Figure 6. Marginalised ID errors from Euclid-like data on each cosmological parameter, for a wCDM model (left), and with the addition of w a (centre) and 
/nl (right), for the local configuration. For each parameter, we show the result obtained using all four observables and their combinations. The error bars in 
red include the Planck priors. We also overlay the current constraints from WMAP7, baryon acoustic oscillations and Hq measurements in light green, and the 
size of the current uncertainty on /nl from the CMB bispectrum in cyan. The results on the standard parameters remain unaffected by the addition of /nl- 
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lensing 

2D clust. phot. 
2D clust. spec. 
3D clust. + shot 
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.11 (.020) 
.079 (.0083) 


.012 (.0016) 
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.017 (.0090) 


.0040 (.00088) 
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.0052 (.0039) 
.0058 (.0049) 


.033 (.0058) 
.044 (.0071) 


.013 (.0033) 
.015 (.0035) 


.011 (.0066) 
.011 (.0076) 


.050 (.041) 
.057 (.051) 


.24(18) 
.28 (.21) 


12 (6.5) 
11 (7.0) 



Table 4. Expected marginalised errors from Euclid-like data for all parameters for the equilateral configuration. The bias is kept scale-independent here. The 
constraints on /nl from galaxy clustering are strongly weakened, as the bias is now scale-independent. The results from lensing remain instead strong. 
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o-(wo) 


cr(w„) 


o"(/nl) 
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2D clust. phot. 
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.32 (.027) 
.34 (.042) 


.13 (.0042) 
.10 (.0041) 
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.82 (.44) 


2.6 (.96) 
3.3 (2.0) 
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12(11) 


lens. + 2D phot. 
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.0082 (.0074) 


.12 (.010) 
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.0094 (.0086) 
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8.6(8.2) 


simple sum 
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.13 (.022) 


.040 (.0041) 


.023 (.018) 


.28 (.22) 


1.6 (.77) 


11(10) 



Table 5. Expected marginalised errors for all parameters for DES for the local PNG case. The numbers within parentheses include the forecasted priors from 
Planck CMB temperature power spectrum. 
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lensing 
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2D clustering, photo-z 
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2D clustering, photo-z 


3.8 (3.6) 


3.2 (3.0) 


2.8 (2.6) 


2D clustering, spectroscopic 
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2D clustering, spectroscopic 


8.6 (7.9) 


6.7 (6.2) 


5.8(5.6) 


3D clust + shot noise 


.13 (.031) 


.12 (.031) 


3D clust + shot noise 


5.7 (5.5) 


4.4 (4.2) 


3.6 (3.5) 


lensing + 2D photo 


.037 (.035) 


.035 (.033) 


lensing + 2D photo 


4.2 (3.7) 


2.8 (2.6) 


2.4 (2.3) 
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.052 (.047) 


.040 (.039) 


lensing + 2D spectro 


6.3 (5.7) 


4.9 (4.6) 


4.1 (3.9) 



Table 8. Forecasted results for non-Gaussianity as a function of the fidu- Table 10. Marginalised forecasted results from Euclid for non-Gaussianity 

cial value of Wq, without and with Planck priors. These results show that as a function of the fiducial value of the bias b^iz), without and with Planck 

a Euclid-like survey should be able to distinguish between a cosmological priors, 
constant and a different dark energy model. 
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2D clustering, photo-z 
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6.7 (6.2) 


8.1 (7.3) 


3D clust + shot noise 
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4.9 (4.6) 
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Table 9. Forecasted results for non-Gaussianity as a function of the fiducial 
value of /nl without and with Planck priors: detection versus upper limit 
scenarios. These results show that a Euclid-like survey will be able to either 
detect models with primordial non-Gaussianity down to the level of /nl ~ 
10, or to find strict constraints around a Gaussian model. 



combined case (lensing + 2D photometric galaxy clustering), the 
uncertainty on / NL degrades from <x(/ N l) = 2.6 to 0"(/nl) = 4.1 
when Planck priors are included. We can see the results sum- 
marised in Table [9] from which we can conclude that this model 
would be detected with high significance by Euclid. 

Galaxy biasing 

The scale-dependent non-Gaussian correction to the galaxy linear 
bias is proportional to the Lagrangian bias b g — \. Therefore, in the 
limit where the scale-independent part is b g (z) = 1, the effect of 
PNG will vanish making it impossible to constrain / N l from galaxy 
clustering on large scales. On the other hand, it is interesting to no- 



tice that in the presence of highly biased tracers, such as luminous 
red galaxies (LRG) or quasars, the effect - and the constraints - are 
maximised. 

We look at how much the constraints are modified if we alter 
our fiducial bias by a fraction /;,, i.e. we take 

b^(z)= i+f b (Vm-1 ) , (47) 

and we consider two cases with f b = i± 0.25. We can see in Ta- 
ble [10] that the constraints on / NL are indeed strongly dependent 
on this choice, as expected, except obviously for the lensing case 
where there is no biasing. 



6.5 Scale-dependent non-Gaussianity 

The local shape of non-Gaussianity can be obtained e.g. when mul- 
tiple scalar fields give a contribution to the curvature perturbations. 
However, in many models of inflation, non-Gaussianity is gener- 
ated in a scale-dependent way. For this reason, a new parameter 
has been introduced by Chen ( 2005 1 in the equilateral case and by 
|Byrnes et al.| ( |2010| > in the local case: the spectral index of non- 
Gaussianity, rif NL , defined so that the effective /nl(&) is written as 

- / k Y" NL 

/nl(*) = /nl , (48) 

where /nl is the value of /nl at the pivot scale A: p i v . We choose the 
pivot scale to be near the logarithmic centre of the expected data 
l |Cortes et al.|2007) and fix it to k p „ = 0.02 fr/Mpc. 
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1.6 (.94) 


37 (19) 
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lens. + 2D phot. 


.014 (.0069) 


.0035 (.00077) 


.0034 (.0032) 


.033 (.0050) 


.011 (.0030) 


.0035 (.0032) 


.037 (.035) 


.19 (.13) 


2.8 (2.6) 


simple sum 


.024 (.011) 
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.0079 (.0066) 


.038 (.0089) 
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.0075 (.0065) 
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.39 (.25) 


2.8 (2.7) 



Table 6. Comparison of summing two matrices obtained for Euclid versus combining them appropriately including all the cross-signals. Results between 
parentheses include Planck priors. 



Euclid data ajfrnJ ^("/nl) 



lensing 


68 (58) 


.66 (.59) 


2D clustering, photo-z 


14 (9.6) 


.38 (.26) 


2D clustering, spectroscopic 


23 (14) 


.64 (.38) 


3D clust + shot noise 


10 (7.6) 


.28 (.21) 


lensing + 2D photo 


5.3 (4.3) 


.18 (.14) 


lensing + 2D spectro 


6.6 (5.2) 


.17 (.12) 



Table 11. Forecasted results from Euclid-like data for scale-dependent non- 
Gaussianity and degraded constraints on the scale-independent part ^jl- 
The fiducial model has here ftjL = 30. 



For the simplest case of local non-Gaussianity, in which a sin- 
gle field is responsible for generating the curvature perturbation 
such as the curvaton scenario (Huang 2010, Byrnes et al. 2011 1, it 
is possible to take the scale dependence of / NL out of the integrals, 
and simply apply it to the bias variation Ab(k) ( Shander a et al.|20lTj 
|Desjacques et al.|201 1| . In this case, by adding to our param- 
eter array, and fixing the fiducial fy L = 30, we obtain the results of 
Table 1 1 1| The best constraint comes from combining lensing, clus- 
tering, and the Planck prior, giving <t(/j/ nl ) ^ 0.12. The constraints 
on /nl get in this case weakened. These constraints are similar to 
to the forecasted limits from the Planck CMB bispectrum ( |Sefusatti] 
|et al.|2009| >, and fully independent. 

6.6 Smooth bias parametrisation 

The constraints obtained so far are very conservative in the treat- 
ment of the bias. By assigning to every redshift bin an independent 
nuisance bias parameter, we are actually allowing for much more 
freedom than it is physically reasonable to expect. In particular, we 
have good reasons to expect the galaxy bias to be a smooth function 
of redshift. Therefore, we propose an alternative, less conservative, 
parametrisation, in which the bias is assumed to be a polynomial 
function of redshift. In the following, we will use a third-order 
polynomial: 

b(z) = b + b x {z - 1) + b 2 (z - If + b 3 (z - l) 3 , (49) 

with fiducial values b = V2,2>i = (2^2)-\b 2 = -(16V2)-\fci = 
(64 VI)- 1 . The expansion in (z - 1) is required if we want the bias 
to approximate Vl + z over the interval < z < 2. Notice that this 
expansion converges only in this interval, and thus it would have to 
be modified for a broader redshift range. We have also checked that 
by expanding up to the third order we are able to approximate the 
Vl + z function to better than 1 % in this range. 

When including the shot noise marginalisation in this case, 
we also parametrise it as a third-order polynomial instead of leav- 
ing it as a completely free set of independent nuisance parameters 
for each bin, for the same reasons of requiring that physical quan- 
tities should have a smooth evolution in redshift. The forecasted 
errors shrink in this case however only marginally, as we can see in 
Table [H 
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.015 
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.0048 


.0047 


.0046 


.0045 


lensing + 2D spectro 


.0066 


.0076 


.012 


.024 



Table 13. Forecasted results from Euclid for the measurement of the galac- 
tic bias. A third-order polynomial expansion is assumed. 



Finally, we evaluated how well could galaxy bias be mea- 
sured with these observables. Referring to the smooth polynomial 
model of Eq. |49j, we calculated the marginalised errors on its 
coefficients, which can be seen for the different Euclid probes in 
Table [T3] Again, the combination of clustering+lensing gives the 
strictest results, since in this case all the other parameters are much 
better constrained, yielding an accuracy on the bias at the level of 
1%. 

6.7 Dependence on the cutoff scales 

Small scales 

Here we study how the constraints change if we vary the maximum 
mode k m . iX , / max used in the forecasts. Since the constraints on each 
individual parameter will not in general vary monotonically, we 
will study the most interesting quantity, the figure of merit (FoM) 
T of the survey with respect to a set of N parameters i?, defined 
similarly to Albrecht & Bernstein ( 2007 ) as 

n*) = — N 1 , (so) 
n,=i o\vi) 

where each <t(#,) is the width of the error ellipsoid along the axis 
defined by the i-th eigenvector of the Fisher matrix. In other words, 
it is easy to obtain their products by taking the determinant of the 
Fisher matrix as 

fW = ViFWI • (51) 

We then study the total FoM, the FoM obtained for the dark energy 
parameters only, and for the dark energy + / NL parameters. We can 
see the result in Fig. [7] where we report the evolution of the FoM 
for lensing and the 3D power spectrum as a function of / max and 
fe max respectively. In the 3D case, we evolve the scale cutoff with 
redshift as explained in the previous sections, by imposing that the 
variance of the density field stays constant in all redshift bins. In 
addition, to avoid considering exceedingly small scales, we impose 
a further condition that in every bin we will only consider scales 
k < = 2k max (z = 0). 

As expected, the results are strongly dependent on the choice 
of the minimum scales used. For the case of galaxy clustering, we 
can see that we could still largely improve the constraints by ex- 
tending the analysis deeper into the non-linear regime, with the use 
of more advanced biasing models reliable on smaller scales. On the 
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.0043 (.00091) 
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.019 (.0030) 
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.089 (.013) 
.026 (.0056) 
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.031 (.014) 
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.27 (.22) 
.085 (.052) 
.13 (.031) 
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.67 (.28) 
.67 (.10) 


3.0 (3.0) 
7.0 (6.6) 
4.4(4.1) 


lens. + 2D phot, 
lens. + 2D spec. 


.014 (.0068) 
.016 (.0071) 


.0035 (.00075) 
.0049 (.00081) 


.0033 (.0031) 
.0038 (.0034) 


.033 (.0049) 
.043 (.0053) 


.011 (.0030) 
.014 (.0035) 


.0034 (.0032) 
.0041 (.0038) 


.036 (.035) 
.039 (.035) 


.18 (.13) 
.21 (.13) 


2.8 (2.6) 

4.9 (4.6) 



Table 12. Expected marginalised errors for Euclid-like data for all parameters, using a parametric form of the galaxy bias. The numbers within parentheses 
include the forecasted priors from Planck CMB temperature power spectrum. For 3D, includes marginalisation over FoG. 
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Figure 7. Dependence of the Euclid-like constraints on the small-scale cuts. 
In the top panels we show the evolution of the full figure of merit as a func- 
tion of the minimum scale considered in the analysis: Z max for the 2D part 
(left) and i max at z. = for the 3D part (right). In the following rows, we 
show the corresponding FoM relative to different combinations of few se- 
lected parameters: wo, w a , /nl- The dashed lines denote Euclid only, while 
the solid lines include the Planck CMB priors (temperature power spectrum 
only). The triangles mark the scale which was chosen in the main analysis. 



other hand we see that for weak lensing we have already included 
most of the signal with our choice of Z max . It can be seen that more 
conservative choices would still yield a very similar result; for ex- 
ample by choosing more conservatively Z milx = 5000, we have that 
the marginalised error on / NL with Planck priors degrade from 19 
to 22 only, and the total FoM degrades by only a factor of 3. 
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Figure 8. Dependence of the Euclid constraints on the large-scale cuts for 
the 2D and 3D cases. The FoM relative to the /nl parameter only is shown. 
The value chosen in the main analysis, k m \ n = 10~ 3 /i/Mpc is marked with 
a triangle. The FoM for the remaining cosmological parameters is largely 
independent from this choice. 

Large scales 

We then also look at the effect of changing the large-scale cutoffs 
hmn, Imm- Due to the larger error bars, this is largely unimportant 
for most cosmological parameters, with the important exception of 
/nl: due to the form of the scale-dependent bias, the large scales are 
the region where this parameter is most constrained, and thus the 
forecasted FoM on /nl is strongly dependent on this choice, as we 
can see in Fig.[8]for both the 2D and 3D cases. However, due to the 
limited size of the surveys, we decide that a conservative choice can 
hardly push to scales larger than k m \ n = 10~ 3 /z/Mpc and / m ;„ = 5, 
which we therefore adopt. A further reason to discard larger scales 
is to avoid the regime where general relativistic corrections become 
important, as described e.g. by Yoo ( 2010| >. 

6.8 Effect of the halo model inaccuracies 

As we have seen in the Section [3] the accuracy of the halo model 
at small scales is only in the range of 10%. For this reason we have 
tested how severely this impacts on the Fisher matrix forecasts de- 
rived using this model. For this purpose, we have calculated two 
sets of Fisher matrices for the Euclid spectroscopic survey, using 
the 3D power spectrum observable. For the first matrix of each set, 
we have derived the power spectra used the halo model; for the sec- 
ond, we have used the approach by Smith et al. (2007). We can see 
in Fig.|9]the results using scales up to /c max = 0.15 (as used in the 
main results) and 0.3 at z = respectively, from which we can con- 
clude that the our results are rather robust in the transition to the 
non-linear regime. 

For weak lensing, where smaller scales are considered, the dis- 
crepancy between the two models for the non-linear matter power 
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Figure 9. Effect of the uncertainties in the modelling of the non-linear 
regime onto the Fisher matrix forecasts for the 3D power spectrum in the 
local configuration. The ellipses in blue and red denote here the forecasted 
posteriors obtained using the halo model and the power spectra from the 
method by Smith et al. 1 2007 1 respectively. Only scales up to /: ma x — 

0.15 

ft/Mpc and i max = 0.3 /i/Mpc at z = were used in the panels. 



spectrum grows larger. Fig.[T]shows that the model by Smi th et aT] 
( |2007| > deviates from N-body simulations by more than 20% for 
k > 2/i/Mpc, whereas the halo model is significantly more accu- 
rate on these scales. Consistently, previous studies have shown that 
|Smith et aT1j2007f underestimates the convergence power spectrum 
by more than 20 % at / > 1000, while the halo model is reliable up 
to / ~ 50,000 (see Fig. 9 in |Hilbert et^pOOg} ). For this reason 
we only consider the halo model in our weak lensing analysis. 

6.9 Redshift binning 

A further question which may be asked is how important is the 
choice of the redshift tomography used, i.e. how much would the 
results change if we used a different binning. We can see in Fig.|10| 
that this choice is not critical for the lensing and the 3D power spec- 
trum, as in these cases the results have largely converged when we 
take more than a few bins. On the other hand, for the projected 2D 
spectra this choice is very important: due to the increasing number 
of cross-correlations between the bins, the signal-to-noise increases 
significantly up to a rather high number of bins. For computational 
reasons, we have decided to carry our analysis using 12 redshift 
bins, where most cases are quite close to saturation. 

6.10 Euclid Red Book 

The detailed specifications of the Euclid mission are still evolv- 
ing. While the current study was performed based on the so-called 
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Figure 10. Dependence of the Euclid constraints on the redshift tomogra- 
phy. We show the evolution of the figure of merit as a function of the num- 
ber of redshift bins considered in the analysis, for the all parameters (top), 
for the dark energy parameters only (middle) and for /ml (bottom panel). 
All these results include the Planck priors, and the choice used in the main 
analysis (M = 12) is marked by a triangle. We can see that the results for 
lensing and the 3D power spectrum have already converged once more than 
a few bins are used, while for the projected 2D cases the increase in number 
of bins brings in much additional signal. 



Yellow Book ( |Laureijs et al.||2009|, the specifi cations have since 
evolved to the Red Book ( |Laureijs et al.|20TT) . In this section we 
describe how our results change as a consequence of this update. 

We show in Table|14|the parameters which have changed from 
the Yellow to the Red book; for the latter, in many cases two val- 
ues are given: minimum requirement and goal. As it can be seen, 
the 'goal' specification is generally the same as the Yellow Book 
value, while the 'requirement' is less ambitious. The main excep- 
tion is the galaxy density in the spectroscopic case, which has in- 
creased due to a deeper flux limit. Table [T3] reports the changes in 
the forecasted constraints for the Red book 'requirement' specifi- 
cations for the local PNG case. The Red Book 'goal' specifications 
are very similar to the Yellow book results presented above. We can 
see that in general the constraints degrade for the photometric sur- 
vey, and slightly improve for the spectroscopic part, mainly due to 
the increased number density in this latter case. 

In addition, we also show the results obtained using a fiducial 
bias derived from semianalytic models of galaxy formation (Orsi 
|et al.|20 10), which are available for the spectroscopic case only. In 
this case the constraints on /nl degrade further, as this tabulated 
bias is approximately 15% lower than the model b(z) = Vl + z 
which we use in our main results. 
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Red Book Update Euclid photometric Euclid spectroscopic 



Parameter 


description 


Required 


Goal 


Required 


Goal 


Mz)/(1 + z) 


redshift uncertainty 


0.05 


0.03 


0.001 


0.001 


z 


median redshift 


0.8 


1.0 


1.0 


1.1 


n 


galaxy density 


30 arcmin~ 2 


40 arcmin~ 2 


1 .20 arcmin~ 2 


1.20 arcmin~ 2 


A 


surveyed area 


15 000 sq deg 


20 000 sq deg 


15 000 sq deg 


20 000 sq deg 


dN/dz(z) 


galaxy distribution 


Smail et al. 


Geach et al. 



Table 14. Changes in the specifications for the Euclid Red Book (Laureijs et al. 201 1 1. The parameters which are not shown are unchanged from Table^ 



Red Book 




o-(fifc) 


o-(n m ) 


a{K) 


cr(n) 


o"(a- 8 ) 


cr(w ) 


o-(w a ) 


o"(/nl) 


lensing 


.078 (.017) 


.019 (.0019) 


.016 (.011) 


.16 (.014) 


.063 (.0041) 


.030 (.015) 


.17 (.12) 


1.0 (.47) 


73 (27) 


2D clust. phot. 


.095 (.036) 


.0087 (.0040) 


.050 (.024) 


.15 (.030) 


.043 (.0039) 


.034 (.026) 


.50 (.31) 


1.9(1.2) 


5.8 (5.5) 


2D clust. spec. 


.076 (.015) 


.0065 (.0021) 


.015 (.011) 


.099 (.015) 


.050 (.0038) 


.031 (.015) 


.090 (.058) 


.54 (.28) 


6.4 (5.9) 


— , tab bias 


.082 (.016) 


.0069 (.0023) 


.017 (.012) 


.11 (.016) 


.056 (.0039) 


.036 (.017) 


.10 (.063) 


.59 (.30) 


8.8 (8.0) 


3D + shot noise 


.067 (.0083) 


.0048 (.0010) 


.023 (.0037) 


.027 (.0064) 


.033 (.0035) 


.042 (.017) 


.13 (.029) 


.59 (.092) 


4.1 (4.0) 


— , tab bias 


.071 (.0084) 


.0051 (.0010) 


.025 (.0038) 


.029 (.0065) 


.036 (.0035) 


.048 (.020) 


.14 (.031) 


.63 (.098) 


6.1 (5.8) 


lens. + 2D phot, 
lens. + 2D spec. 

— , tab bias 


.030 (.0087) 
.030 (.0077) 
.031 (.0080) 


.0056 (.00097) 
.0052 (.00089) 
.0054 (.00092) 


.0048 (.0044) 
.0051 (.0037) 
.0054 (.0040) 


.050 (.0065) 
.048 (.0058) 
.051 (.0061) 


.016 (.0036) 
.016 (.0036) 
.018 (.0037) 


.0052 (.0048) 
.0054 (.0045) 
.0060 (.0051) 


.054 (.052) 
.051 (.035) 
.055 (.039) 


.32 (.20) 
.35 (.15) 
.38 (.16) 


4.7 (4.5) 
5.7 (5.3) 
7.9(7.3) 



Table 15. Expected marginalised errors from Euclid for the Red Book specifications (minimum requirements) for the local non-Gaussianity case. The numbers 
within parentheses include the forecasted priors from Planck CMB temperature power spectrum. The 3D clustering includes marginalisation over the shot 
noise parameters in each redshift bin. For the spectros copic part of the s urvey, we show both the results obtained using our standard form for fiducial bias 
b(z) °c Vl + z and the results from the tabulated bias by Orsi et al. 2010 (in italic). 



7 CONCLUSIONS 

In this work we have studied to what accuracy will the two-point 
statistics of future surveys of the large-scale structure determine 
the cosmological parameters, with a particular focus on the non- 
Gaussianity parameter /nl, considering the most relevant local, 
equilateral and orthogonal bispectrum configurations. We have per- 
formed a Fisher matrix analysis using the specifications of the up- 
coming Dark Energy Survey and a Euclid-like survey based on the 
Euclid Assessment Study Report, using both its spectroscopic and 
photometric parts. We have chosen these surveys as examples of 
the future DETF stages III and IV, and the results are likely compa- 
rable with other surveys in the same class, such as in particular the 
planned American mission WFIRST. 

We have combined all the relevant data sets, including their 
covariances. In particular, we have considered the projected 2D 
galaxy power spectrum, always including curvature and using the 
exact calculation, discarding the Limber approximation on large 
scales, as this introduces inaccuracies which are important in the 
presence of scale-dependent galaxy bias from primordial non- 
Gaussianity. 

We obtained that the strictest constraints on /nl are expected 
from the combination of weak lensing and photometric galaxy clus- 
tering; in this case we find for the local case cr(/ NL ) =^ 3 for Euclid 
and - 8 for DES, when also including priors from the tempera- 
ture power spectra of the Planck CMB mission. In the cases of or- 
thogonal and equilateral configurations, the constraints from galaxy 
clustering degrade greatly, due to the reduced scale-dependence of 
the galaxy bias, while the constraints from weak lensing remain 
at a similar level. Finally, the constraint on scale-dependent non- 
Gaussianity is for the local case <r(nf NL ) = 0.12 when the fiducial 
scale-independent part is f NL = 30. The level of these constraints is 
comparable to the expectations from the Planck CMB bispectrum, 
and fully independent (Komatsu & Spergel 2001; Sefusatti et al. 
|2009l >. 

We have also studied the effect of updating the description of 



Euclid to the latest Red Book specifications; in this case using the 
Red Book 'goal' parameters leaves the forecasts largely unchanged, 
while we found that using the 'requirement' parameters degrades 
the constraints on PNG to cr(/ NL ) ^ 5. 

Further extensions of this approach, which will be useful not 
only for forecasting but also for the likelihood analysis of the up- 
coming real data, will include on one hand expanding the total 
covariance of the two-point statistics to include observations from 
clusters of galaxies and the correlations with the cosmic microwave 
background; on the other hand, this approach will be extended to 
the inclusion of the three- and four-point statistics, which will be 
instrumental in the search for the higher-order gNL and tnl infla- 
tionary parameters. 
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APPENDIX A: PRIMORDIAL NON-GAUSSIANITY AND 
THE LARGE-SCALE STRUCTURE 

We review here the effects of primordial non-Gaussianity on the 
large-scale structure, focussing on the results which are used in this 
analysis. 
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Al Bispectra and skewness 

The bispectrum of the Bardeen potential O in the local, equilateral 
and orthogonal cases is given by < |Taruya et al.| 2008, Sch midt &| 
|Kamionkowski|20 1 0) 

B^(k u k 2 ,k 3 ) - 2/ NL + 2perms.] , (Al) 



xF(ki , At) F(k 2 , At) F(k, At) , 
The skewness is finally defined for all cases as 
S 3 (M) = (6 3 M )/(6 2 M ) 2 = (6 3 M )/<r 4 . 



(A8) 



(A9) 



flT(*i.*2,fe) - 6/= q L {-[P*(fc 1 )P*(fc 2 ) + 2perms.] 

-2 [Po(h)Po(k 2 )P 9 (.k 3 )f 3 

+ [pV\h)P%\k 2 )Po(k3) + 5 perms.] } ,(A2) 

Bg(k u k 2 ,k 3 ) - 6/*{-3[/'o(^)/ > ofe) + 2perms.] 

-8[Po(^i)/ > *(fe)/'*fe)] 2/3 

+ 3 [p; /3 (fe)4 /3 (fe)P*(fe) + 5perms.]}(A3) 

It is interesting to compute the lowest moments of the correspond- 
ing linear density field 6: the variance and the skewness. The vari- 
ance is defined as 



<r 2 (M) = (S 2 M ) 



2n 2 J 



dkk 2 P(k)F 2 (k, At), 



(A4) 



where we have introduced the smoothed field 6 M (k) = 6(k)F(k, At), 
and F(k,M) is a filter function of mass resolution At. We use a 
top-hat function in real space. 

The third momentum of the smoothed density field can be 
written in terms of the three-point function 



C d 3 ki r d 3 k, r d 3 k 3 
J Wj (2^<^^^(k 2 )^(k 3 )>. (A5) 



<<&> = 



By substituting the bispectra from Eqs. l |Al||A2||A3[ l, and integrat- 
ing the Dirac delta to give k3 = -ki - k 2 so that we can define 
k = |ki + k 2 |, this expression can be simplified for the three cases 
extending the calculation by Desjacque s et al.|(2009| > 

^ fdk^ r^ffe) 

2n 2 ) 2 J 'aft) Jo 2 a( 



(sir 



(2n 



*(k 2 ) 



I' 



x d/ia(k) 



1+2 



P(k)a 2 (k 2 ) 



~P(k 2 )a 2 {k) 
xF(k { , At) F(k 2 , At) F(k, At) , 



(A6) 



(2 



4gj f dhk\a{h) f dk 2 kla(k 2 ) 
■ n ) Jo Jo 

I dfia(k) 



aHhWikj 



+ 2 perms. 



-2 



P{h)P(k 2 )P(k) 



2/3 



tfik^aP-ihWQc) 
/ >1/3 (yt 1 )P 2/3 (fc 2 )P(/fc) 



a 2 ' 3 (ki)a 4 ' 3 (k 2 )a 2 (k) 
xF(k u M) F(k 2 , M) F(k, At) , 



+ 5 perms. 



(A7) 



/ c3 \ ort 



3fi 
(2. 



* 2 ) 2 Jo 
X J djia(k) 



dk\k\a(k\) I dk 2 k\a{k 2 ) 
Jo 

' P(h)P(k 2 ) 



P{h)P{k 2 )P{k) 



aHk^aHh) 

2/3 



+ 2 perms. 



+3 



a 2 (/ti)a 2 (/t 2 )a 2 (^) 



a 2 l 3 (h)tfl 3 {k 2 )a 2 (k) 



+ 5 perms. 



A2 Mass function 

The halo mass function n(M, z) gives the number density of haloes 
of mass M at redshift z. Various analytical methods have been de- 
veloped to compute this quantity starting from the statistical prop- 
erties of Gaussian primordial perturbations, starting with Press & 
|Schech ter ( 1 974 1 . Fitting formulae that improve agreement with N- 



body simulations have also been presented ( Sheth & Tormen|1999 



[20021 |JenkTns et al.|200l]|Warren et al.|2006UTinker et al.|2008 



|Pillepich et al.|2010| >. 

The halo mass function is modified by the presence of PNG 
([Matarrese et al.||200"0{ |Lo Verde et al.||2058{ |Maggiore & Riotto] 
2010a, [Lam & Sheth|2009[ l. Since analytical models are based on 
a set of simplistic assumptions, they most robustly predict the ratio 
between the mass function generated by non-Gaussian and Gaus- 
sian initial conditions with the same power spectrum (see |Giannan-| 
|tonio & Por ciani ( 2010J for a comparison of the different methods 
against N-body simulations). 

For this reason, in our forecasts we use the fitting formula 
given by PPH08 for the Gaussian mass function and multiply it 
by the factor \Lo Verde et al.|2008| > 



^LV f^i/NL^ 



1 + 



S 3 <T 2 (q 



r4 



2qS 2 c 



6 y[qS c 
dS 3 (cr) a I gS] 
dhur 2 \ cr 2 



-1 + 



(A10) 



in the non-Gaussian case. Here S c - 1.686 is the threshold for the 
collapse of a linear density perturbation and the corrective factor q, 
of order unity, is a heuristical correction which may be applied to 
the collapse threshold to improve the agreement with N-body sim- 
ulations, and which we set to unity for the reasons given in Section 
[2] Note that PNG enters the expression for the halo mass function 
through the linear density skewness S 3 (M) which is linearly pro- 
portional to /nl. 



A3 Halo bias 

The halo linear bias factors may be obtained from the expression 
for the halo mass function using the peak-background split formal- 



ism (Bardeen et al. 


1986||Cole & Kaiser| 1989| |Mo & White|1996[ 


|Catelanetal.|1998 


i. In the Gaussian case, the bias is obtained by 



taking the logarithmic derivative of the mass function with respect 
to the collapse threshold. This gives a mass and redshift dependent 
bias coefficient. 

Recent generalisations of the peak-background split to cases 
with PNG have shown that the corrections to the mass function 
generate two extra terms in the expression for the halo bias: one of 
them is a small addition to the Gaussian bias while the second one 
introduces a scale-dependent term (see Eqs. |7} and (9} in the main 
text) (Dalai et al.|2008l|Slosar et al.|2008[patarrese & Verde|2008[ 
Afshordi & Tolley 2008; Giannantonio & Porciani 2010, Schmidt 
& Kamionkowski 2010; Desjacq ues et al.|2011} . 
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APPENDIX B: PLANCK FISHER MATRIX 

The Planck Fisher matrix was kindly made available to us by 
Jochen Weller. It was calculated according to|Albrecht et al.| ( |2009| l; 
|Rassat et al. 1(2008) and used in |Sartoris et al.|(2010| >. As"described 
in Section 3.8 of Sartoris et al. ( 2010) , it is computed adopting the 
parameterisation i? = (a> m ,& s ,lnAs, £J&,n s ,T): with A s the ampli- 
tude of primordial perturbations, the size of the sound horizon 
at the last-scattering surface, and r the optical depth to Compton 
scattering which has been marginalised over. The calculation of the 
Fisher matrix is based on the conservative assumption that only the 
143 Ghz channel is used for science analysis. In this case, the beam 
size is i?fwhm = 7.1', and the temperature and polarisation sensitiv- 
ities are a T = 2.2/jK/K and o> = 4.2^/K/K, respectively. In order 
to avoid polarisation foregrounds, only multipoles with / min = 30 
have been considered together with a sky fraction f s ^ y = 0.8 to min- 
imise the effects of galactic foregrounds. The transformation of the 
Fisher matrix for the parameter set (S1„,,Q. A , h, <x 8 , flj,, Wo, w a , n s ) is 
done as in Rassat et al. (2008). It is important to highlight that the 
Planck Fisher matrix does not include any information on primor- 
dial non-Gaussianity, as CMB constraints on PNG can be calcu- 
lated from higher-order statistics only (especially the bispectrum). 
So adding these priors will not improve the constraints on / NL . 
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